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HASTI 

A  Numerical  Calculation  of  Two-Dimensional 
Lagrangian  Hydrodynamics  Utilizing  the  Concept 
of  Space-Dependent  Time  Steps 


by 

Phillip  L.  Browne  and  Martha  S.  Hoyt 


ABSTRACT 


A  numerical  experiment  is  described  in  which  an  attempt  was  made  to 
reduce  the  calculation  time  of  a  two-dimensional,  Lagrangian,  hydrody¬ 
namics  code,  without  loss  of  accuracy.  The  reduction  was  accomplished 
by  utilizing  a  space-dependent  time  step;  i.e.,  the  frequency  with  which 
a  point  is  calculated  depends  on  the  rate  at  which  changes  are  occurring 
in  the  vicinity  of  the  point.  This  system  permitted  savings  to  be  made 
by  eliminating  calculations  in  regions  undergoing  relatively  slow  changes. 
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I.  INTRODUCTION 


In  many  numerical  models  approximating  physical  or  mathematical 
systems  where  successive  approximations  are  used  to  describe  model 
changes  as  the  calculation  proceeds,  there  are  often  regions  in  which 
changes  are  taking  place  more  rapidly  than  in  others.  The  usual  com¬ 
puting  procedure  is  to  extrapolate  all  regions  of  the  system  at  the 
interval  required  by  the  region  of  most  rapid  change,  resulting  in  un¬ 
necessary  calculations  for  regions  undergoing  slow  changes.  It  was  the 
purpose  of  this  program  to  see  if  eliminating  these  superfluous  calcula¬ 
tions  might  save  calculation  time,  with  little  loss  in  accuracy.  The 
disadvantages  inherent  in  such  a  procedure  are  the  additional  storage 
required  and  the  increased  complexity  of  the  code  logic,  neither  of 
which  turns  out  to  be  serious  in  practice,  especially  when  a  modem 
computer  with  a  large  memory  is  used. 

Because  numerical  calculation  of  two-dimensional  hydrodynamics  is 
time-consuming  and  relatively  complicated,  and  because  the  authors  are 
somewhat  familiar  with  two-dimensional,  Lagrangian  hydrodynamics,  this 
field  was  chosen  for  experimentation.  It  was  decided  to  experiment 
with  a  Magee-type  code  inasmuch  as  it  had  been  studied  in  some  detail 
by  one  of  the  authors  and  his  coworkers.1  The  Magee  code  was  initiated 

p 

at  Los  Alamos  by  H.  Kolsky,  and  has  been  developed,  improved,  and  used 

3 

extensively  by  S.  R.  Orr  and  others.  It  is  appropriate  to  mention 

here  that  some  other  work  using  a  space-dependent  time  interval  has  been 

4 

done  by  Goad,  using  other  methods  in  the  code  WAT  with  good  success. 

The  concept  of  a  space -dependent  time  step  can  easily  be  applied 
to  one -dimensional  hydrodynamics  problems;  and  the  authors  feel,  with 
some  conviction,  that  the  same  approach  can  be  used  with  problems  in- 
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volving  heat  or  radiation  diffusion,  in  either  one  or  two  dimensions. 

It  is  possible  to  speculate  further  and  suggest  that  similar  schemes 
be  used  in  numerical  systems  wherever  successive  approximations  are 
used. 

II.  A  BRIEF  DESCRIPTION  OF  THE  MAGEE  CODE  (NON-SPACE -DEPENDENT  TIME 

INTERVAL) 

Consider  a  physical  system  of  cylindrical  symmetry,  which  is  repre¬ 
sented  in  a  cylindrical  (r,z)  coordinate  system  by  a  mesh  consisting  of 
quadrilateral  zones  (as  shown  in  Fig.  l),  rotated  through  an  angle  of 
one  radian  about  the  z  axis.  (The  use  of  this  angle  simplifies  the 
formulae  and  calculations  by  eliminating  the  factor  2rr.)  Thus,  each 
quadrilateral  zone  represents  an  element  of  volume.  Because  of  the 
cylindrical  symmetry,  all  quantities  are  dependent  only  on  r  and  z. 
Inasmuch  as  the  Lagrangian  viewpoint  has  been  adopted,  the  mesh  is  con¬ 
sidered  to  be  imbedded  in  and  to  move  with  the  fluid. 


Fig.  1.  A  Typical  Magee  Mesh 


In  order  to  describe  a  hydrodynamic  system  mathematically,  it  is 
necessary  to  use  the  principles  of  conservation  of  mass,  momentum,  and 
energy,  together  with  an  equation  of  state  relating  the  thermodynamic 
quantities.  Also,  it  is  necessary  to  know  the  initial  condition  of  the 
system  and  some  boundary  conditions  as  the  hydrodynamic  activity  unfolds. 
In  the  following  paragraphs,  the  required  variables  and  difference- 
approximations  of  these  equations  will  be  enumerated. 

IIA.  The  Variables 

In  Magee,  many  of  the  quantities  are  assumed  to  be  associated  with 
the  points  (the  intersections  of  grid  lines)  and  are  identified  by  a  pair 
of  integers,  (i,j),  as  subscripts,  i.e.: 

\ 

The  coordinates  of  a  point,  (r.  .  and  zJ  .) 

i  ,y 

The  velocity  components  of  a  point,  (r.  .  and  z  )  )  (l) 

The  acceleration  components  of  a  point,  (r  .  and  z.  .) 

1,0  i>0 

Other  quantities  are  assumed  to  be  associated  with  the  zones  and  are 
identified  by  a  pair  of  fractional  subscripts,  (i-^,j-|),  i.e.: 


The  pressure  and  Richtmyer-Von  Neumann  artificial  dissipative \ 

5  6 

term,  *  (p+q),  i  .  x  =  PlJL 

The  internal  energy  per  unit  original  volume,  ^(2) 

The  mass,  density,  and  relative  volume,  m.  x  .  ,  , 


The  stability  number,  (w/b ) 


i-ibo-ir 


The  times,  t,  t+At,  t+At/2,  at  which  quantities  are  as  sinned  to  be 
known  or  calculated  are  identified  by  superscripts  such  as  n,  n+1  and 
n+§.  The  superscript  n+|  is  also  used  to  identify  the  change  in  a 
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quantity  between  tn,  tn+"*".  The  superscript  0  denotes  the  initial  value 
of  a  quantity  at  the  start  of  the  problem. 

IIB.  The  Difference  Equations 

The  conservation  of  mass  is  not  expressed  explicitly  in  a  separate 
equation;  but,  rather,  it  is  assumed  implicitly  in  writing  the  momentum 
and  energy  equations.  Because  the  internal  energy  equation  is  applied 
to  a  zone  as  an  entity,  for  that  equation  the  mass  of  each  zone  is 
assumed  to  be  conserved.  And,  because  the  momentum  equation  is  applied 
to  a  point  as  an  entity,  for  that  equation  the  masses  associated  with 
each  point  are  considered  to  be  always  the  same  (specifically,  one-fourth 
of  the  mass  of  each  of  the  adjacent  zones). 

To  simplify  this  discussion,  the  difference  equations  will  be  listed 
in  a  plausible  sequence  for  calculation,  not  necessarily  the  sequence 
actually  used  in  Magee  (in  which  a  number  of  variations  have  been  tried) . 
Several  passes  through  the  mesh  might  be  required  to  perform  the  calcu¬ 
lations  as  presented  here. 

Assume  that  all  calculations  have  been  performed  up  through  cycle  n 
(t=tn),  so  that  for  all  points  and  zones  there  are  available  the  values 
of: 


n 


n 


.n-h 

r  ei. 


.n- 


r  v11 

>  »  } 


Pn, 


and  e 


n 


(3) 


1  3 

The  accelerations  at  time,  n,  are  calculated  from  an  approximation  ’ 
of  the  law  of  conservation  of  momentum  given  in  Eq.  (4).  The  arguments 
used  to  derive  these  equations  are  too  lengthy  to  be  given  here,  but  will 
hopefully  be  derived  in  a  subsequent  report.  The  quantities  used  in  Eq. 
(4)  are  associated  with  a  group  of  points  and  zones  as  shown  in  Fig.  2. 
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These  accelerations  my  now  be  used  to  advance  the  velocities  of 
all  points  to  time,  n+|,  by 

r"*  =  rn^  +  rn[i(tntl+tn)  -  KtV1)] 

jn+t  .  pi  +  gn[i(tnrt+tn)  .  (5) 


Then,  these  velocities  are  used  to  advance  the  coordinates  of  all  points 
to  time,  n+1,  by 
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The  next  steps  involve  the  calculation  of  V,  P ,  € ,  etc . ,  for  the 
zones  at  time,  n+1. 

The  actual  volume  \sQ  of  a  quadrilateral  zone  can  be  calculated 
exactly  as  the  sum  of  the  volumes  of  two  zones  of  triangular  cross 
section  (see  Fig.  3). 

r 

f 


z 

Fig.  3*  The  Two  Triangular  Zones  Making  up  a  Quadrilateral  Zone 
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(10) 


TIn+l  .  n+1  /.  0 
Vd  =  1D  Ad 

the  I®  having  been  obtained  at  the  start  of  the  problem  in  order  to 
calculate  the  masses  of  the  zones  by 


0,0 

Pn 


(ID 


These  masses,  assumed  to  be  constant  in  time,  have  been  used  in  Eq.  (4). 
If,  at  any  time,  the  density  of  a  zone  is  required,  it  is  obtained  by 


.  n  m 
n+1 _ n_ 

P°  “  ,  n+1 


0 


m. 


(12) 


5,6  . 

vr  7  -T  t 


The  Richtmyer-Von  Neumann  artificial  dissipative  term *  is  cal- 
culated  from  the  rate  of  volume  change  for  each  zone  by 


where 


n+1 


/,  0%2  .n+1  1  f  AV3^ 

p0  (vn«}2  V  / 


Av 


n+i; 


=  vn+1  -  v11  <  0, 


and 


(13) 


,  n+1 


area  of  quadrilateral  at  time,  n+1. 


In  the  case  of  a  zone  expansion,  i.e.,  >  0, 


n+1 

q  =0 


(13a) 


The  pressure,  p,  and  internal  energy,  e,  of  each  zone  are  obtained 
from  the  equation  of  state 

pn+1  -  p(v"V+:l)  (14) 
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and  from  the  equation  of  conservation  of  energy,  d£  =  -PdV,  in  difference 
form. 


en+1  =  en  -  |<P%n+V+1)AVB+;'  (15) 

If  (l4)  and  (15)  cannot  he  solved  simultaneously  for  the  two  unknowns, 
n~bl  n+1 

p  and  €  ,  then  iterative  techniques  are  used. 

7 

The  White  stability  number  for  each  zone  is  calculated  by 


i+t 


r  r  w  'N2-in+1  1  +  3pn+V1+1  ^Atn+K2  j 

l  vlX'1  j  - - ^+1  ■'vTTS/  +  4  TH+I 

p0  A  1  T 


(16) 


n+3 /2 

and  the  time  interval.  At  '  ,  to  be  used  for  the  next  cycle  is  chosen 
so  that  the  maximum  of  the  stability  numbers  for  all  the  zones  with 


satisfies 


or 


Atn+3^  -  At"* 


S' 

0.035  <  j—  <  0.14 


(IT) 


max 


0.22^5  <  v  <  0.^5 

max 

Frequent  calculation  of  the  total  energies  in  regions  of  a  problem 
are  useful,  although  not  necessary.  The  internal  energy  is  added  by  zones 


<“>"  ■  L  ^-M4  ^4,34 
1,3 


(18) 


and  the  kinetic  energy  is  added  by  points 
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III.  THE  HASTI  CODE:  GENERAL  DESCRIPTION 

Because  the  HAST I  code  was  experimental,  emphasis  was  placed  on 
achieving  a  code  that  would  work  properly  and  accurately,  yet  be  con¬ 
venient  to  change.  Less  attention  was  given  to  conserving  storage  and 
machine  time  by  the  use  of  short  cuts  and  clever  manipulation.  Conse¬ 
quently,  there  may  be  instances  where  the  code  is  inefficient  in  these 
respects . 

In  any  problem,  the  points  of  the  mesh  can  be  classified  into 
regions,  or  groups,  according  to  the  maximum  time  intervals  at  which 
the  quantities  in  those  regions  may  be  calculated  without  producing  in- 
stability  or  inaccuracy.  Goad  has  used  a  method  of  listing  points  in 
a  set  of  tables.  We  have  chosen  a  different  method,  in  which  two  time 
values  are  stored  at  the  end  of  a  cycle  for  each  point: 

last  time  at  which  the  point  was  calculated 
next  time  at  which  the  point  should  be  calculated  (20) 

There  is  also  available  a  time 

=  minimum  of  all  t^+J  (21) 

which  is  the  next  time  a  calculation  must  be  made  for  some  points  in 
the  problem.  Furthermore,  there  is  stored  a  time 
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tn+'1'  =  next  time  at  which  all  points  must  be  calculated 

major 


(22) 


The  purpose  of  this  time  is  to  bring  the  calculation  together  for  analysis 
(pointing,  plotting,  etc . ) ,  and  its  choice  is  arbitrary.  Variations  in 
time  intervals  were  made  by  powers  of  two. 

For  each  point,  i, j,  among  the  quantities  stored  are:* 


vn+l 

n 

a.  . 


representing  r.  . ,  z.  at  t.  . 
representing  r^,  at 

representing  r.  ,,  z  .  at  t“ 

J-J  «J  tJ  ^9  J 


>(23) 


^4,14  represelItine  vi4>3-4’  *i4.H  at  *i,4 

At+  -  .  i  representing  the  next  time  interval  allowed  by 
the  stability  number 


To  illustrate  the  method,  consider  the  problem  to  be  at  the  stage 
where,  in  the  logical  mesh,  there  is  one  region  of  points  (denoted  • 
in  Fig.  4)  for  which  tn+1  =  tj^.  All  other  points  have  tn+1  >  tjjjj. 

The  values  of  tn  for  points  may  differ.  The  authors  believe  that  the 
arguments  and  methods  to  be  presented  also  apply  to  cases  where  there 
are  many  such  regions  and/or  regions  of  more  conplex  shapes.  And,  in¬ 
deed,  the  problems  run  with  the  code  have  turned  up  no  exceptions  to 
this  theory. 

Because  the  •  points  are  probably  in  a  region  undergoing  more  rapid 
changes  than  the  surrounding  points,  this  rapidly  changing  domain  is  en¬ 
closed  by  a  "buffer  region",  one  zone  thick  (denoted  by  the  shaded  area 
in  Fig.  4) .  Enough  calculations  are  made  for  the  buffer  region  so  that 
the  stability  number  can  be  calculated  for  these  zones  in  order  to  alert 
the  code  if  the  disturbance  is  moving  into  the  buffer  region,  ^buffer 
zone  in  Fig.  4  is  delineated  by  three  kinds  of  x  points  (\  i, «) 

-"■For  sinplification  of  notation  in  the  subsequent  discussion  and  flow 
diagrams,  several  quantities  will  be  represented  by  the  same  symbol. 

IT 


Fig.  4.  A  Single  HA.STI  Logical  Mesh,  Shoving  the  Various  Types 
of  •  ,  x,  and  o  Points 

and  by  five  kinds  of  o  points  (  o  ,  p,  o-»  ,0,  c/)  all  of  which  are 
defined  by  a  succession  of  tests  on  t  ,  as  shown  in  the  flow  diagram 
(Fig.  5). 

In  order  to  have  the  required  quantities  available  at  the  proper 
time,  the  code  makes  two  passes  through  the  mesh,  taking  points  in  the 
order 
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Fig.  5*  HA.STI  -  First  Pass-Through  Mesh 
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(  START  ^ 


Fig.  6.  HA.STI  -  Second  Pass-Through  Mesh 
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i  =  1,  j  -  1,  2,  .  .  .  Jlagt 
i  =  2,  j  =  1,  2,  .  .  .  3 last 


1  =  J  =  1;  2>  •  •  •  J^ast 

On  the  first  pass  (see  Fig.  5)  calculations  are  made  for 

rn+l,  Vn+1,  wn+1,  and  At+ 

while  on  the  second  pass  (see  Fig.  6),  the  quantities  determined  are 

an+\  yn+3/2 ^  revj_se(3_  tn,  and  tn+^  (26) 

Because  of  space  limitations,  the  flow  diagrams  do  not  include  the 
equations  used  for  calculations ;  these  equations  are  part  of  the  dis¬ 
cussion  of  the  problem  illustrated  in  Fig. 

IV.  THE  FIRST  PASS 

As  the  calculation  proceeds  on  the  first  pass  through  the  mesh  of 
the  example  in  Fig.  4  [following  the  ordering  scheme  shown  in  Eq.  (2*0  ], 
circle  (o)  points,  dot  (•)  points,  and  x  points  will  be  encountered. 

The  first  pass  is  used  to  calculate  the  coordinates  and  thermodynamic 
quantities . 

IVA. 

In  the  exampl e  in  Fig.  4,  the  first  points  in  the  mesh  requiring 

calculation  are  the  circle  (o)  points.  The  coordinates  of  these  points 

must  be  advanced  to  ta.  because  this  information  is  required  to  obtain 

min 

accurate  volumes,  V,  for  the  adjacent  dot  (•)  points.  Failure  to  ad¬ 
vance  the  circle  (o)  points  can  and  has  resulted  in  the  undesirable 
effect  of  an  associated  zone  apparently  compressing  when,  actually,  it 


21 


is  expanding,  thereby  producing  an  artificial  viscosity  by  Eq.  (13)  when 
there  should  be  none.  Such  a  fictitiously  produced  perturbation  has  a 
tendency  to  spread  through  the  mesh  in  the  form  of  oscillations  of 
points  about  their  proper  positions,  and  may  result  in  instability  if 
the  stability  condition  is  not  made  more  severe. 

Consider  the  pictorial  time  scale  of  Pig.  7,  in  which  the  times  for 
a  dot  point  are  represented  below  the  line  by  tn  and  t^J^j  and  the  times 
for  an  adjacent  circle  point  are  represented  above  the  line  by  (t /)n  and 
(t')n+1. 


(t') 


n 


(t") 


n-f-1/2 


(t') 


n+l/2 


(t') 


n+1 


' 

o 

• 

t 


n 


t 


n+1 

min 


Pig.  7.  A  Pictorial  Time  Scale  for  Adjacent  Circle  and  Dot  Points 


The  coordinates,  rn,  and  accelerations,  a11,  of  the  circle  point  are 

known  at  ( t ,)nj  the  velocities,  vn+s,  are  known  at  (t7)11*^  =  ^[(t ')n 

+  (t7)  ]:  and  it  is  desired  to  advance  the  coordinates  to  ta.  .  The 

n+i  niin 

simplest  way  to  do  this  is  to  assume  that  v  *  is  valid  for  the  whole 
interval  from  (t')n  to  (t7)n+1  and,  hence,  by  Eq.  (6): 


rO  - 


n 


n+irj_n+l 


(27) 
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n+l 

If  more  accuracy  is  desired,  the  velocity  can  be  adjusted  from  (t ')  2 

to  (t")n4g  =  |[(t')n+t^3  by  an  equation  similar  to  Eq.  (5): 

(■ v")n+ *  -  vn+^  +  an[(t")n+^-(t')n+2)  (28) 


before  using  Eq.  (27). 
n*f*l 

Once  r(t  .  )  has  been  used  to  calculate  the  related  volumes,  it  is 
nun 

then  deadvanced  to  r  by  Eq.  (27);  and  the  velocity  is  corrected  back  to 

n4vl 

v  2  by  Eq.  (28),  so  that  the  proper  values  are  in  storage  for  the  two 
times  (t')n  and  (t')n  1  still  associated  with  that  point.  In  order  to 
delay  the  corrections  until  after  the  volumes  have  been  calculated,  the 
corrections  are  made  after  other  x  or  dot  points,  further  along  in  the 
mesh,  have  been  calculated  (see  Fig.  5)*  The  delay  is  necessary  in 

nj  -i 

order  to  have  the  coordinates  of  all  corners  of  a  zone  at  t  .  ,  so  that 

min 

the  volume  can  be  calculated  more  accurately. 

IVB.  The  Dot  (»)  Points 

Calculations  for  the  dot  (•)  points  are  quite  straightforward. 
Knowing  tn,  t^4^,  a11,  and  v11,  one  uses  Eq.  (6)  to  compute 


rft”!1)  = 
v  min7 


n  ,  n+4/^n+l 

r  +  v 


(29) 


and,  subsequently,  to  compute  V(t^4^) ,  w(t^T^) ,  etc.  by  Eqs.  (10),  (13), 
(l4),  (15)  and  (16).  The  ftt+  is  the  largest  At  that  satisfies  Eq.  (17). 
This  procedure  is  equivalent  to  what  is  normally  described  as  time¬ 
splitting  and  -doubling. 

IVC .  The  x  Points 

The  x  points  are  buffer  points  similar  to  the  circle  points;  but 
whereas  the  circle  points  are  deadvanced  to  their  original  positions 
after  use,  the  x  points  are  left  in  their  new  positions. 

Consider  the  pictorial  time  scale  in  Fig.  8  in  which  times  for  a 
dot  point  are  represented  above  the  line  by  tn  and  t**^;  and  the  times 
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for  an  adjacent  x  point  are  represented  below  the  line  by  (t')n  and 

<t')n+1. 


t 


n 


n+1 


t 
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/  /Nn 

(t) 


(t7 


n+1/2 


(t'> 


n+]/2 


(t"') 


n+1/2 


(t') 


n+1 


Fig.  8.  A  Pictorial  Time  Scale  for  Adjacent  x  and  •  Points 

The  coordinates,  rn,  and  accelerations ,  a11,  of  the  x  point  are  known  at 
(t')n.  The  velocities,  vn+^,  are  knovn  at  (t')11*5  *  ?j[(t ,)n+('t')Il+13i 
and  it  is  desired  to  advance  the  coordinates  to  t“1y).  Again,  as  a  first 
approximation  one  can  use  Eq.  (27)  or,  for  a  better  approximation,  Eqs. 
(28)  and  (27).  Now,  however,  it  is  desired  to  retain  r(t".^) ,  so  the 
velocity  must  be  advanced  to  (t'")n+^  =  ^[t^^+(t/)n+1]  by 

(V'0n+^  =  (v")n+^  +  an[(t///)n+^-(t/')I1+^]  (30) 


V.  THE  SECOND  PASS 

At  the  start  of  the  second  pass,  the  coordinates  and  thermodynamic 

el+1 

quantities  have  been  calculated  for  all  the  necessary  points  at  . 
During  this  pass  it  is  necessary  to  calculate  accelerations  and  velocities, 
and  to  analyze  the  stability  number  in  order  to  set  up  a  new  pair  of  times, 
where  necessary. 
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New  accelerations  at  t^  are  calculated  only  for  the  dot  points 
by  Eq.  (4) .  Before  new  velocities  are  confuted,  the  times  must  be  ana¬ 
lyzed. 

For  each  point,  the  At"*"  *s  of  the  four  adjacent  zones  are  examined; 

+ 

and  the  smallest,  nonzero  At  is  chosen  to  obtain 


-  <*E>  + 


Ln+lv 


(31) 


which  is  the  next  time  at  which  that  point  must  be  calculated.  The 
point  is  thus  controlled  by  the  adjacent  zone  in  which  the  most  severe 
stability  condition  is  present#  The  (t.  . )  value  is  then  compared 

X^J  101x1  n+P 

to  8X1(1  the  smaller  of  the  two  is  chosen  as  the  actual  t“  .  to 

be  used  for  the  point. 

For  the  dot  points,  the  new  velocities  are  then  calculated  as  in 

Eq.  (5), 


Tn+3/2  .  Ta^~  +  .  i^n+l^n  }] 

v  min'  2V  i,j  min'  2V  min  i,j' 


=  vn+£  +  aCt^1)  [|(tn+2-tn  .)] 
v  min'  2'  i,j  i,j'J 


(32) 


The  remainder  of  Fig.  6  deals  with  the  shifting  of  times  in  storage  to 
prepare  for  the  next  cycle. 


VI.  RESULTS 

There  is  no  need  to  describe  in  detail  the  problems  used  to  test  the 
code.  Suffice  it  to  say  that  the  problems  included  a  number  of  one-dimen¬ 
sional  systems,  such  as  planes,  cylinders,  and  spheres,  in  order  to  in¬ 
sure  that  HASTI  did  not  affect  the  symmetry  of  such  systems.  Further, 
several  simple  problems  with  real  two-dimensional  motions  were  con?>leted. 

The  problems  involved  a  variety  of  situations  involving  shocks,  rarefactions, 
free  surface  motion,  sliding  motion  along  a  fixed  boundary,  and  fixed 
boundaries . 
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In  all  cases ,  tiae  coordinates  and  densities  of  the  problems  run  by 
the  HASTI  code  agreed  with  the  problems  run  by  a  Magee -like  code  to 
within  one  percent.  In  accurate  tests  of  machine  time  used  (excluding 
input-output),  the  HASTI  problems  ran  two  to  three  times  faster  than  the 
Magee  ones.  Because  only  a  limited  effort  was  spent  making  HASTI  an 
efficient  code,  while  the  Magee-type  code  is  very  efficient,  there  is  a 
good  possibility  of  increasing  this  time-saving  factor. 

VII.  CONCLUSIONS 

The  time-saving  possibilities  of  a  space -dependent  time  step  in  the 
numerical  calculation  of  two-dimensional,  Lagrangian  hydrodynamics  have 
been  illustrated;  and  a  fairly  straightforward  method  has  been  suggested. 

The  method  apparently  does  not  save  as  much  time  as  the  system  used  by 

4 

Goad,  but  it  is  felt  that  further  work  and  experimentation  with  HASTI 
could  yield  greater  time  savings.  If  less  emphasis  were  placed  on  accu- 
racy,  (for  example,  eliminating  some  of  the  velocity  adjustments  in 
Section  IV)  and  if  more  emphasis  were  placed  on  efficient  use  of  the 
machine  (saving  quantities  calculated  on  one  point  for  use  with  another 
point),  farther  time  savings  could  result.  The  method  was  found  to  be 
sufficiently  general  and  flexible  so  that  changes  and  refinements  of  the 
code  could  be  made  without  too  much  difficulty.  For  exanple,  a  modified 
version  of  the  Magee  ” fudge"  was  added  to  stop  any  point  that  attempted 
to  cross  the  diagonal  of  the  zone. 
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